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Abstract. Energy in nuclear matter is, in practice, completely characterized at different densities 
and asymmetries, when the density dependencies of symmetry energy and of energy of symmetric 
matter are specified. The density dependence of the symmetry energy at subnormal densities pro- 
duces mass dependence of nuclear symmetry coefficient and, thus, can be constrained by that latter 
dependence. We deduce values of the mass dependent symmetry coefficients, by using excitation 
energies to isobaric analog states. The coefficient systematic, for intermediate and high masses, is 
well described in terms of the symmetry coefficient values of a% — (31.5— 33.5) MeV for the vol- 
ume coefficient and a s a = (9-12) MeV for the surface coefficient. These two further correspond 
to the parameter values describing density dependence of symmetry energy, of L ~ 95 MeV and 
tf sym ~25MeV. 
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INTRODUCTION 

Symmetry energy is first encountered in nuclear physics in the empirical nuclear energy 
formula as the term describing drop in nuclear binding with increasing neutron-proton 
asymmetry. That term is quadratic in the asymmetry, expressing the charge symmetry, 
symmetry of nuclear interactions under neutron-proton interchange. In nuclear matter, 
the energy as a function of neutron and proton densities may be expanded in the matter 
asymmetry. Knowledge of density dependence of the expansion coefficient and of the 
energy of symmetric matter suffices in practice for the determination of nuclear energy 
and pressure at any density and asymmetry, which is essential for calculations of neutron 
star structure [1]. The coefficients of expansion for finite nuclei and nuclear matter are 
connected in a nontrivial manner. 

Specifically, the density dependence of the symmetry-energy expansion coefficient in 
nuclear matter gets tied to a variation with mass of the symmetry-energy coefficient in 
the empirical energy formula, due to changing contributions to the energy of nuclear vol- 
ume and of surface where nuclear density drops. As a consequence, it may be possible to 
constrain the density-dependence of symmetry energy using the mass-dependence of the 
symmetry coefficient. In the following, we deduce the symmetry energy coefficient as a 
function of nuclear mass, using nuclear energies. Moreover, we employ Skyrme-Hartree- 
Fock (SHF) calculations for half-infinite matter and for spherical nuclei, combined with 
Hohenberg-Kohn functional theory, to constrain the density-dependence of symmetry 
energy. 



SYMMETRY COEFFICIENT AND ISOVECTOR DENSITY 



Considerations 10, d] of the competition between nuclear surface and volume in storing 
neutron-proton asymmetry leads to the following formula for nuclear energy 

E (N, Z) = -a v A + a s A 2 ' 3 + (N-Z) 2 + a c -^+ E mic , (1) 

where E m i c represents microscopic contributions to the energy and where the mass- 
dependent symmetry coefficient a a (A) follows from 

A A A 2 / 3 



a a {A) a v a a s a 

In the above, and a s a are the volume and surface symmetry coefficients and Eq. © 
states that the capacity of the system for asymmetry is a sum of the interior and surface 
capacitances. 

The energy per nucleon in nuclear matter, expanded in the relative asymmetry r\ = 
(p n -Pp)/p, becomes 

j = ^(p) + ^(p,r 1 )^(p)+S(p)7 1 2 - O) 
The coefficient S is typically expanded around the normal density as 



S(p)=a v a +—(p-p ) + -^(p-po) 2 + .... (4) 
JpO lopQ 

Considerations [4] of the nuclear Hohenberg-Kohn (HK) energy functional |0], 
smoothed out to make it analytic, imply that the mass -dependent symmetry coefficient 
may be generally represented as 



A _ 1 
a a {A) a v a 

where the isovector density p a is defined as 



drp fl (r) , (5) 



P.(r) = * (f '- P ' (,) . (6) 

Here, r\y is asymmetry in the interior of matter. 

In the limit of weak Coulomb effects, the isovector and isoscalar densities change only 
within the second order in asymmetry across an isobaric chain [0]. That implies that, in 
the continuum limit across medium and heavy nuclei, the proton and neutron densities 
can be in practice described in terms of just 4 parameters, radius parameter, difusenesses 
for isoscalar and isovector densities and the relative displacement of isovector and 
isoscalar densities. The last 2 parameters are tied to the density dependence of symmetry 
energy. 
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FIGURE 1. Isoscalar (solid curves) and isovector (dashed) density in half-infinite nuclear matter, 
together with a local approximation to the isovector density (short-dashed), as a function of position, 
for sample Skyrme interactions. 



In the limit of short nonlocality range in the symmetry part of the HK energy func- 
tional, compared to nuclear densities, the isovector density is 

In Fig. [0 we show densities from our SHF calculations of half-infinite nuclear mat- 
ter. Before such calculations have been carried out Farine et al. [6] . For p > p /4, the 
exact isovector densities exhibit Friedel oscillations around the respective local approx- 
imations given by Eq. ©. This implies that differences between proton and neutron 
densities, and even proton densities alone, can directly express the density dependence 
of S at p > p/4. At p < p/4, nonlocalities in the HK energy functionals become strong 
as evidenced in the discrepancies between exact and approximate isovector densities. 

From the energy of half-infinite matter and/or vector density we can determine as 
and a s a coefficients in the energy formula (OQ) for the Skyrme interactions. While the 
surface coefficients, ag ~ 19MeV, vary little between different interactions, the surface 
symmetry coefficients a s a vary widely, from 9 to 75 MeV. The coefficients a s a are tightly 
correlated to the symmetry-energy slope parameter L, see Fig. [21 A correlation also is 
found between scaled L and K iym . Those correlations are likely due to the lack of other 
energy scale in the energy functional than pr> 
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FIGURE 2. Correlation between the slope parameter L and inverse of the surface symmetry parame- 
ter a s a , for Skyrme interactions from the compilation by Stone [7]. The parameters are scaled with a^. 



a a {A) FROM ISOBARIC ANALOG STATES 

Competition between different physics terms within an energy formula hampers the 
ability to learn about mass dependence of the nuclear symmetry coefficient by fitting 
directly a formula to the ground-state nuclear energies yL |8fl. However, the unwanted 
competition may be practically eliminated by generalizing the energy formula to lowest 
states of given net isospin within a nucleus y90, amounting to the replacement of the 
symmetry term in (QQ): 

a a (A) [ — ] = 4a a (A) -j- Aa a (A) [ -± ] . (8) 

This replacement absorbs the so-called Wigner term from Z?nuc- With the formula gen- 
eralization, it becomes possible to deduce the symmetry coefficient within a single 
nucleus, by using excitation energies to the isobaric analog states (IAS) representing 
ground states of neighboring nuclei [80, with 

AE = 4a a (A) ^ + A£ mic . (9) 

A 

We use the tabulated energies of isobaric analog states 111011 and microscopic correc- 
tions to energies by Koura et al. to deduce symmetry coefficients for individual nu- 
clear masses. Those corrections include deformation effects. Deduced coefficient values 
range from a a ~ 9MeV for light A < 10 nuclei to a a ~ 22.5MeV for A > 200. Figure[3] 
shows inverse coefficient values plotted against inverse cube root of mass number. For 
A > 20, the value systematic is approximately linear, as expected from Eq. ©. A fit 
with Eq. ©, produces coefficient values of = 32.9 MeV and a s a = 11.3MeV. Similar 
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FIGURE 3. Inverse of asymmetry coefficient, in the left scale, as a function of inverse cube root of mass 
number, in the bottom scale. For convenience, top and right scale show, additionally, the mass number and 
the coefficient before inversion. The filled circles represent coefficients extracted from IAS excitation 
energies with the microscopic corrections applied. The line shows a fit to the IAS data for A > 20. 
The crosses and squares represent, respectively, the coefficients from spherical Sk5 SHF calculations 
and coefficients from a Thomas-Fermi model that best describes the data. 



coefficient values are obtained when trying to describe the a a (A) data in terms of the 
Thomas-Fermi theory Q3L [8J] or spherical SHF |fl2l1 . 

The parameters arrived at in different ways are further summarized in Fig. |4l The 
Thomas-Fermi results might be less reliable than other, because model yields an evo- 
lution of density profile in the surface region with mass, which is not fully sup- 
ported by electron scattering data. Overall, the IAS data suggest parameter values of 
a% = (31.5 - 33.5) MeV for the volume coefficient and a s a = {9- 12) MeV. Using the 
correlations between ci^ and a s a coefficient values and L and £sym> sucn as in Fig. [2l 
within the Skyrme model, we further arrive at values of L ~ 95 MeV and K sym ~ 25 MeV. 
These represent a nearly linear S(p) at p < po- 

Further efforts in this investigation strive to understand better the Coulomb effects 
for finite nuclei [8] and to exploit skin sizes as well as details in p p (r) for further 
constraining 5(p) at p < po. 
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FIGURE 4. Symmetry-energy parameter values in the plane of a v a ja s a and a v u , from comparing theoret- 
ical expectations to data. 
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